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The production of ^®A1 in massive stars is sensitive to the ^^Na(Q;,p)^®Mg cross section. Recent 
experimental data suggest the currently recommended cross sections are underestimated by a factor 
of ~40. We present here differential cross sections for the ^^Na(a,p)^®Mg reaction measured in the 
energy range Ec.m. = 1.7 — 2.5 MeV. Concurrent measurements of Rutherford scattering provide 
absolute normalisations which are independent of variations in target properties. Angular distribu¬ 
tions were measured for both po and pi permitting the determination of total cross sections. The 
results show no significant deviation from the statistical model calculations upon which the recom¬ 
mended rates are based. We therefore retain the previous recommendation without the increase in 
cross section and resulting stellar reaction rates of a factor of 40, impacting on the ^®A1 yield from 
massive stars by more than a factor of three. 


The observation of ^®A1 in the galactic medium, 
through 7 -ray emission from its daughter nucleus ^®Mg, 
provided direct evidence for ongoing nucleosynthesis in 
the galaxy [1]. While the origins of ^®A1 remain the sub¬ 
ject of discussion, the C/Ne convective shell within mas¬ 
sive stars is a candidate site [2]. A sensitivity study of the 
reactions influencing A1 production in massive stars has 
indicated a significant dependence on the ^^Na(a,p)^®Mg 
reaction rate, which acts as a proton source for the 
^®Mg(p, 7 )^®Al reaction [3]. Specifically, it was found that 
an increase in the ^^Na(a,p)^®Mg rate by a factor of 10 
would lead to an increase in ^®A1 production by a factor 
of 3. 

The ^^Na(a,p)^®Mg rate adopted in Ref. [3] is obtained 
from statistical model calculations. While earlier exper¬ 
imental data do exist [4, 5], these were excluded due to 
a lack of understanding of the target properties during 
the intense beam bombardment. As a consequence there 
are significant uncertainties in the experimentally deter¬ 
mined resonance strengths. 

A recent direct measurement of the reaction cross sec¬ 
tion in inverse kinematics was made to resolve these 
experimental uncertainties [6]. A ^^Na beam was inci¬ 
dent on a gas cell containing ^He and outgoing protons 
corresponding to the ground and first-excited states in 
^®Mg detected. The cross sections measured in the re¬ 
gion Ec.m. = 1-7-2.5 MeV were ~40 times greater than 
statistical model calculations. Such an increase is signif¬ 
icantly larger than that required to alter the production 
of ^®A1 by a factor of 3. 

A similar, although less dramatic, disagreement with 
statistical model calculations has been reported for the 
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^^S(q;,p)^®C 1 reaction [7]. It is noted in Ref. [8] that 
the measured cross sections significantly exceed the ex¬ 
pected single particle strength and that, in light of the 
^^Na(Q;,p)^®Mg results also, there is an urgent need for 
additional {a,p) data in the 20 < A < 50 region. 

In this letter we report on a new measurement of the 
^^Na(Q;,p) cross section in forward kinematics covering 
the energy range Ec.m. = 1-7-2.5 MeV. Our methodology 
exploits the simultaneous detection of Rutherford scat¬ 
tered a particles to remove dependencies on properties 
of the target, such as thickness and stoichiometry, which 
has impacted previous measurements. Discussions of this 
methodology may be found in, for example. Refs. [9, 10]. 

Measurements were made at the Aarhus University 5- 
MV Van de Graaff accelerator. A schematic of the ex¬ 
perimental setup is shown in Fig. 1. A ^He beam with 
laboratory energy between 1.99 and 2.94 MeV was used 
to bombard a carbon-backed NaCl target. The beam 
was stopped 70-cm downstream of the target position in 
a suppressed Faraday cup connected to a current inte¬ 
grator. Typical beam currents were in the range 200- 
500 ppA. 

Two double-sided silicon strip detectors (DSSDs) were 
mounted in the scattering chamber to provide energy and 
angle information for outgoing charged particles. A 322- 
/im annular DSSD was mounted upstream of the target, 
covering laboratory angles between 140° and 163° and a 
40-/rm-thick, quadratic DSSD provided coverage at labo¬ 
ratory angles between 60° and 120°. The annular detec¬ 
tor was mounted with the junction side, which has a 4/im 
dead layer, facing the target. In this orientation the dead 
layer acts as a degrader foil, increasing the energy sep¬ 
aration between backscattered a particles and protons 
populating ^®Mg. Protons populating the ground and 
first excited states in ^®Mg were sufficiently energetic to 
punch through the quadratic DSSD and were stopped. 
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FIG. 1. A schematic of the experimental setup within the 
scattering chamber. The incoming "^He beam is indicated by 
the arrow. The NaCl target was orientated at 45° relative to 
the beam axis. Two double-sided silicon strip detectors were 
used to detect outgoing protons and a particles, see the text 
for details. For clarity both front- and back-side segmentation 
of the detectors is shown. 


and unambiguously identified, in a 1500-/im silicon pad 
detector. 

The target was prepared at Aarhus University by evap¬ 
orating NaCl onto a 10 /ig/cm^ carbon foil. The beam 
energy-loss in the target was calibrated using alpha par¬ 
ticles backscattered from the carbon backing into the an¬ 
nular DSSD. As the target is rotated through 180° the 
energies are shifted due to losses within the NaCl layer 
(see, for example. Ref. [11] for details of this technique). 
A thickness of 65 keV at a beam energy of 3 MeV was 
determined. It should be noted that the target was tilted 
at 45° to the beam axis during all other measurements 
giving an effective thickness of between 92 and 115 keV 
for the range of beam energies used. 

During the experiment elastically scattered alpha 
particles were continuously measured in the quadratic 
DSSD. For pure Rutherford scattering the elastic yield is 
a product of the target thickness and incident beam cur¬ 
rent. This removes any uncertainties due to changes in 
the target thickness or stoichiometry, in addition to un¬ 
certainties in the integration of beam charge. The a scat¬ 
tering data presented in Ref. [12] demonstrate that elastic 
scattering from Na is well described by the Rutherford 
formula for beam energies up to 3 MeV, which covers the 
entire range of measurements here. This is supported by 
a measurement of the angular distribution for elastically 
scattered a particles from Na measured with our setup, 
shown in Fig. 2, which shows excellent agreement with 
Rutherford scattering. 

Rutherford scattering data were also collected for Cl 
throughout the experiment. A comparison between the 
relative amounts of Na and Cl in the target shows no de¬ 
viation from a ratio of 1:1 to within 10% for the duration 
of the experiment, see Fig. 2. Repeat measurements of 
the target thickness were also consistent, indicating no 
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FIG. 2. (a) The measured angular distribution of elasti¬ 

cally scattered alpha particles from ^^Na at a beam energy 
of 2.94 MeV. The solid line is the distribution expected for 
pure Rutherford scattering, (b) The stoichiometric ratio of 
Na and G1 in the target foil, as determined by Rutherford 
scattering, as a function of the integrated beam impinging on 
the foil. 


significant changes in the target properties during the 
experiment. This is not surprising given the relatively 
low beam currents employed, three orders of magnitude 
lower than those used in Ref. [5] where significant target 
degradation was observed. 

Energy spectra for the annular DSSD and quadratic 
DSSD plus pad detector telescope are shown in Fig. 3. 
Proton yields were extracted for both po and pi tran¬ 
sitions across the full energy range covered. Differen¬ 
tial cross sections were obtained using the normalisation 
provided by Rutherford scattering of a particles into the 
quadratic DSSD. 

Examples of measured angular distributions are shown 
in Eig. 4. To permit total cross sections to be determined 
measured differential cross sections were fitted using a 
sum of even-termed Legendre polynomials. This assumes 
a distribution symmetric around dc.m. = 90°) which is 
expected when the cross section is dominated by com¬ 
pound nucleus formation. In the measurements reported 
in Ref. [4] only a single, relatively minor resonance was 
found to exhibit forward-backward asymmetry in the en¬ 
ergy region covered here. Nonetheless, a conservative 
20% uncertainty on the total cross section is assumed. 

At beam energies below ^2.2 MeV, some fraction of 
Pi protons reach the pad detector with insufficient en¬ 
ergy to be registered. The result is a decrease in the 
detection efficiency which is not easily quantified. For 
the two data sets collected below this energy the pi data 
collected using the pad detector are therefore not used 
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FIG. 3. Representative energy spectra from the annular 
DSSD (S3) and the quadratic DSSD plus pad detector tele¬ 
scope (Wl). In the latter case a coincidence between the two 
detectors is required to remove the background due to a par¬ 
ticles stopping in the quadratic DSSD. The effective centre of 
mass energy in each case is given in the plot, see the text for 
details. The small peak at * has an energy consistent with 
the po transition in ^®F(Q;,p)^^Ne, and may therefore be in¬ 
dicative of a thin layer of ^®F on the target surface, see the 
text for further discussion. 


for the fitting of angular distributions. Instead only data 
from the annular DSSD, which suffers no decrease in de¬ 
tection efficiency, is used and an isotropic angular distri¬ 
bution assumed. Applying the same procedure to the pi 
data sets at higher energy results in a decrease in total 
cross section of between 10% and 30%. A 30% uncer¬ 
tainty is assumed for the two lowest pi data points to 
reflect this. 

In two of the eight measurements a weak peak was 
observed ~ 200 keV in energy below po (see the middle 
panel of Fig. 3). This peak may be indicative of a thin 
layer of fluorine on the target surface since it lies at the 
approximate energy expected for ^®F(Q;,po)^^Ne. Under 
these circumstances there may be a contribution from 
^®F(a,pi)^^Ne also, which would not be resolved from 
^^Na(a,pi)^®Mg. The only possible effect of this could 
be to increase the observed cross section. Based upon 
the data presented in Ref. [13] it is estimated that this 
contribution should always be below 10% of the total 
peak yield. Out of caution an additional 10% uncertainty 
is therefore assumed on the lower bound of all pi cross 
sections. 

Protons populating higher lying states in ^®Mg were 
not observed in this work due to the background from 
scattered beam. The contribution from these states to 
the total reaction cross section at the energies measured 



FIG. 4. Proton angular distributions from the ^^Na(a,p)^®Mg 
reaction. The energies given are the effective centre of mass 
energies, corrected for the target thickness. For the lowest 
energy pi data shown reliable differential cross sections could 
only be obtained in the annular DSSD detector, located at 
backward angles in the laboratory frame. The dotted lines 
show fits of Legendre polynomials to the data. 


in this work is expected to be minor due to the reduced 
penetrabilities. In Ref. [6] a Hauser-Feshbach calculation 
is reported which indicates negligible contribution from 
P 2 within the Gamow window (Ec.m. — 1.2-2.2 MeV). 

To account for energy losses within the NaCl layer of 
the target, the measured cross sections are associated 
with an effective beam energy. This is calculated using 
an energy dependence for the cross section as given by 
the statistical model code Non Smoker [14]. The resulting 
effective energies are within 15 keV of the beam energy 
at the target mid-point for all measurements. 

Total cross sections for pg and pi are presented in Ta¬ 
ble I. These values are plotted in Fig. 5 together with 
results from Ref. [6] and the statistical model code Non 
Smoker [14]. We find a significant discrepancy with the 
results reported in Ref. [6], these values being consis¬ 
tently an order of magnitude higher than measured here. 
We can offer no explanation for this discrepancy, how¬ 
ever it cannot be accounted for by the form of the an¬ 
gular distributions assumed in Ref. [6], where data were 
only obtained backwards of Oc.m. = 160°. In the narrow 
angular range between 0c.m. = 165° and 170° where over¬ 
lapping differential cross section measurements exist, the 
absolute values again differ by at least an order of mag¬ 
nitude. It is again worth noting that the absolute nor¬ 
malisation in the present work is provided by Rutherford 
scattered beam from the ^^Na component of the target 
itself. Combined with the relative simplicity of the exper¬ 
imental setup, this provides an extremely robust method 
for the determination of absolute cross sections. 

The Non Smoker results reproduce the measured cross 
sections extremely well in terms of both trend and 
magnitude. The only significant deviation is found at 
Ec.m. = 2.16 MeV and can be understood in terms of 
the strong individual resonance reported in Ref. [5] at 
Ec.m. = 2.14 MeV. If the energy dependence of the Non 
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TABLE I. Angle-integrated cross sections for the po and pi 
branches of the ^®Na(a,p)^®Mg reaction. The final column 
gives the ratio of the measured cross section to that calculated 
using the statistical model code Non Smoker. 


Ec.m. (keV)" 

a Po (mb) 

a Pi (mb) 

u(po +Pi)/aN.s. 

1744 

0.05 (1) 

0.06 it'i) 

1.50 

(29) 

1831 

0.09 (2) 

0.20 (t®) 

2.09 

(46) 

1998 

0.08 (2) 

0.24 (t®) 

0.81 

(13) 

2071 

0.20 (4) 

0.52 itll) 

1.19 

(19) 

2139 

0.28 (6) 

2.42 (ttl) 

3.20 

(58) 

2328 

0.28 (6) 

1.52 (til) 

0.84 

(14) 

2400 

0.57 (11) 

1.59 (1^^) 

0.73 

(11) 

2469 

1.62 (33) 

2.97 (1®°) 

1.18 

(17) 


^ Effective energy corrected for energy loss within the target. See 
the text for details. 



FIG. 5. Cross sections for the reaction ^^Na(Q:,p)^®Mg. The 
energies given are effective energies, corrected for energy 
losses within the target. See the text for details. For compari¬ 
son cross sections from the statistical model code Non Smoker 
[14] and the measurement reported in Ref. [6] are also shown. 

Smoker results is fixed and only the absolute magnitude 
allowed to vary we find a scaling factor of 0.96±0.06 is 
required to best fit our data. 

In conclusion, we have presented cross sections for the 
^^Na(Q!,p)^®Mg reaction in the region E^.m. = 1-74 to 2.47 
MeV. The overall trend and magnitude of the cross sec¬ 
tion are in general found to be very well reproduced by 
the statistical model code Non Smoker. The results are 
also largely consistent with the previous measurements of 
Whitmire et al. [5] and Kuperus et al. [4], though in gen¬ 
eral slightly higher than their results, whereas our mea¬ 
surement is inconsistent with the recent measurement by 
Almaraz-Calderon et al. [6]. 

As mentioned, the only significant discrepancy be¬ 
tween the Non-Smoker statistical model and our mea¬ 
surement is at the energy of the strongest (a,p) reso¬ 
nance at Ec.m. = 2.14 MeV, a resonance which is par¬ 
ticularly strong in the pi channel. From the difference 
between the observed cross sections around 2.07 MeV 
and 2.14 MeV centre of mass energy, we estimate the pi 


and Po resonance strengths for this resonance to be wyi = 
1000(300) eV and wyo = 42(13) eV respectively. Based 
on these resonance strengths, the corresponding single¬ 
resonance contribution to the reaction rate is shown in 
Fig. 6 compared to the Non-Smoker reaction rate. The 
contribution from this resonance in itself exhausts up to 
50% of the Non-Smoker reaction rate (at 2 GK), and 
could therefore potentially increase the total reaction rate 
beyond that of the Non-Smoker rate. At the most im¬ 
portant temperature, 1.4 GK, the temperature at termi¬ 
nation of convective shell G/Ne burning [3], the single¬ 
resonance contribution to the reaction rate is 35% of the 
Non-Smoker reaction rate, with a reduced contribution 
below that temperature. Based on this, we would still 
recommend usage of the Non-Smoker reaction rate for 
^^Na(a,p)^®Mg reaction in astrophysical scenarios, rather 
than the reaction rate indicated in Ref. [6]. The error on 
the reaction rate as evaluated from our experimental data 
is significantly reduced to the level of 30% relative error 
on the reaction rate, except in the temperature region 
around 2 GK where the contribution from the resonance 
could increase the reaction rate by up to 50% as shown in 
Fig. 6, with a corresponding increase in the upper limit 
on the reaction rate. 



FIG. 6. The single-resonance contribution to the total rate 
obtained from Non-Smoker calculations based upon the mea¬ 
sured strength of the resonance at Ec.m. = 2.14 MeV (see text 
for details). 


In summary, we therefore conclude that the reaction 
rate in the key temperature region, around 1.4 GK, 
is consistent with that of the statistical model (Non- 
Smoker), to within approximately 30%. Based on this, 
the resulting ^®A1 production in massive stars as pre¬ 
sented in Ref. [3] still stands. From the results of this 
sensitivity study, in which a 30% ^®A1 production in¬ 
crease is found for a rate-increase of a factor of two, 
the uncertainty in the ^®A1 production corresponding 
to our reaction-rate uncertainty of 30% is expected to 
be at most 10-20%. This level of precision in the 
^^Na(Q!,p)^®Mg reaction rate should therefore be suffi¬ 
cient for detailed comparisons of observed and simulated 
astrophysical ^®A1 production. 
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